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We propose an alternative, nonsingular, cosmic scenario based on gravitationally induced particle 
production. The model is an attempt to evade the coincidence and cosmological constant problems 
of the standard model (ACDM) and also to connect the early and late time accelerating stages of 
the Universe. Our space-time emerges from a pure initial de Sitter stage thereby providing a natural 
solution to the horizon problem. Subsequently, due to an instability provoked by the production of 
massless particles, the Universe evolves smoothly to the standard radiation dominated era thereby 
ending the production of radiation as required by the conformal invariance. Next, the radiation 
becomes subdominant with the Universe entering in the cold dark matter dominated era. Finally, 
the negative pressure associated with the creation of cold dark matter (CCDM model) particles 
accelerates the expansion and drives the Universe to a final de Sitter stage. The late time cosmic 
expansion history of the CCDM model is exactly like in the standard ACDM model, however, there 
is no dark energy. The model evolves between two limiting (early and late time) de Sitter regimes. 
All the stages are also discussed in terms of a scalar field description. This complete scenario is fully 
determined by two extreme energy densities, or equivalently, the associated de Sitter Hubble scales 
connected by pi/pf = (Hi/Hf) 2 ~ 10 122 , a result that has no correlation with the cosmological 
constant problem. We also study the linear growth of matter perturbations at the final accelerating 
stage. It is found that the CCDM growth index can be written as a function of the A growth index, 
YA — 6/11. In this framework, we also compare the observed growth rate of clustering with that 
predicted by the current CCDM model. Performing a x 2 statistical test we show that the CCDM 
model provides growth rates that match sufficiently well with the observed growth rate of structure. 

PACS numbers: 98.80.-k, 95.36.+X, 98.80.Es 



I. INTRODUCTION 

The existence of a cosmological constant A which can 
be used in order to explain the recent cosmic accelera- 
tion, has brought the following major theoretical prob- 
lem: within the framework of the quantum field theory 
the vacuum energy density is more than 120 orders of 
magnitude larger than the observed A value measured by 
the current cosmological data. This is the so called "old" 
cosmological constant problem (CCP) [3, 0. The "new" 
problem Q asks why is the vacuum density so similar 
to the matter density just now? Many solutions to both 
theoretical problems have been proposed in the literature 
• An easy way to overpass the above problems is to 
replace the constant vacuum energy with a dark energy 
(DE) that evolves with time. However the nature of DE 
is far from being understood. Indeed a main caveat of 
this methodology is the fact that the majority of the DE 
models appeared in the literature are plagued with no 
physical basis and/or many free parameters. 

Nevertheless, there are other possibilities to explain 
the present accelerating stage. In particular, the inclu- 
sion of the backreaction in the Einstein Field Equations 
(EFE) via an effective pressure (which is negative for an 
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expanding space-time) opened the way for cosmological 
applications. In these models, the gravitational produc- 
tion of radiation or cold dark matter provides a mecha- 
nism for cosmic acceleration as earlier discussed in 0-Q . 
As a consequence, several interesting features of cosmolo- 
gies where the dark sector is reduced due to the creation 
of CDM matter have been discussed in the last decade 
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In brief, the merits of the particle creation scenario 
with respect to the usual DE ideology are a) the former 
has a strong physical basis namely nonequilibrium ther- 
modynamics, while the latter (DE) has not and b) the 
particle creation mechanism unifies the dark sector (dark 
energy and dark matter), since a single dark component 
(the dark matter) needs to be introduced into the cosmic 
fluid and thus it contains only one free parameter. In- 
terestingly, from the viewpoint of a statistical Bayesian 
analysis models which include only one free parameter 
should be preferred along the hierarchy of cosmological 
models [l6|. We would like to emphasize here that the 
only cosmological models (to our knowledge) which sat- 
isfy the above statistical condition are: 

i) the concordance ACDM which however suffers from 
the coincidence and fine tuning problems 

ii) the braneworld cosmology of (l~7j which however 
does not fit the SNIa+BAO+CMB(shift-parameter) data 
(see Q). 

(iii) the current particle creation model which simul- 
taneously fits the observational data and alleviates the 
coincidence and fine tuning problems [lH, [l5[ . 
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In this paper, we are proposing a new cosmological 
scenario which is complete in the following sense: all the 
accelerating stages of the cosmic evolution are powered 
uniquely by the gravitational creation of massless (at the 
very early stage) and massive cold dark matter particles 
(at the late stages). 

In our scenario, the Universe starts from an unstable de 
Sitter dominated phase (a oc e Hlt ) powered by the pro- 
duction of massless particles, and, as such, there is no the 
horizon problem. Subsequently, it deflates and evolves to 
the standard radiation phase (a oc t 1 / 2 ) thereby ending 
the creation of massless particles. Due to expansion, the 
radiation becomes subdominant with the Universe enter- 
ing in the cold dark matter (CDM) dominated era, in 
which the linear growth of matter fluctuations is taking 
place in a natural way. Finally, the negative pressure 
associated with the creation of cold dark matter parti- 
cles accelerates the expansion and drives the Universe to 
a final de Sitter stage. In addition, the transition from 
Einstein-de Sitter (a oc rj 2 / 3 ) to a de Sitter final stage 
(a oc e Hft ) guarantees the consistence of the model with 
the supernovae type la data and complementary observa- 
tions, including the growth rate of clustering. A transi- 
tion redshift of the order of a few (exactly the same value 
predicted by ACDM) is also obtained. 

The paper is structured as follows. In Sec. II, we dis- 
cuss the basic ideas underlying the particle production 
in an expanding Universe and set up the basic equations 
whose solutions describe the complete evolution of our 
model. In Sec. Ill, we study the linear growth of per- 
turbations, whereas in Sec. IV, we constrain the growth 
index through a statistical analysis involving the latest 
observational results. In Sec. V, a possible scalar field 
description for all stages is discussed, and, finally, in Sec. 
VI we summarize the basic results. 



II. A COMPLETE COSMOLOGICAL 
SCENARIO WITH GRAVITATIONAL PARTICLE 
PRODUCTION 



The microscopic description for gravitationally in- 
duced particle production in an expanding Universe be- 
gan with Schrodinger's [19j seminal paper, which referred 
to it as an alarming phenomenon. In the late 1960s, this 
issue was rediscussed by Parker and others [20T422T ] based 
on the Bogoliubov mode-mixing technique in the context 
of quantum field theory in a curved space-time described 
by the general relativity theory [23| . Physically, one may 
think that the (classical) time varying gravitational field 
works like a 'pump' supplying energy to the quantum 
fields. 

In order to understand the basic approach, let us con- 
sider a real minimally coupled massive scalar field <p 
evolving in a flat expanding Friedman-Robertson- Walker 
(FRW) geometry. In units where K = ks = c = 1, the 
field is described by the following action 



g ap d a 



„»2 /2 

m <p 



(1) 



In terms of the conformal time r\ (dt = 0(77)^77), 
the metric tensor g^ v is conformally equivalent to the 
Minkowski metric r/^, so that the line element is ds 2 = 
a 2 (r])r]^dx IJ, dx 1 ' , where a(rf) is the cosmological scale fac- 

a(i])~ 1 x, one obtains 



tor. Writing the field cf>(r], x) 
from the above action 



X"-V 2 X + mV-- x = 



(2) 



where the prime here denotes derivatives with respect 
to r\. Notice that the field % obeys the same equation of 
motion as a massive scalar field in Minkowski space-time, 
but now with a time dependent effective mass, 



l eff 



(V) 



(3) 



This time varying mass accounts for the interaction be- 
tween the scalar and the gravitational fields. The energy 
of the field x is n °t conserved (its action is explicitly time 
time dependent), and, more important, its quantization 
leads to particle creation at the expense of the classical 
gravitational background [20M22j | . 

On the other hand, in the framework of general rel- 
ativity theory, the scale factor of a FRW type universe 
dominated by radiation (a oc t 1 / 2 ) satisfies the following 
relation [ad + a = 0], or, in the conformal time, a" = 0. 
Therefore, for massless fields (m = 0), there is no par- 
ticle production since Eq. ([2]) reduces to the same of a 
massless field in Minkowski space-time, and, as such, its 
quantization becomes trivial. This is the basis of Parker 
theorem concerning the absence of massless particle pro- 
duction in the early stages of the Universe. Note that 
Parker's result does not forbid the production of mass- 
less particles in a very early de Sitter stage (a" ^ 0). 
Potentially, we also see that massive particles can always 
be produced by a time varying gravitational field. As we 
shall see, such features are incorporated in the scenario 
proposed here. 

In principle, for applications in cosmology, the above 
semiclassical results have three basic difficulties, namely: 

(i) The scalar field was treated as a test field, and, 
therefore, the FRW background is not modified by the 
newly produced particles. 

(ii) The particle production is an irreversible process, 
and, as such, it should be constrained by the second law 
of thermodynamics. 

(iii) There is no a clear prescription of how an irre- 
versible mechanism of quantum origin can be incorpo- 
rated in the EFE. 

Later on, a possible macroscopic solution for these 
problems was put forward by Prigogine and coworkers Q 
using nonequilibrium thermodynamics for open systems, 
and by Calvao, Lima & Waga Q through a covariant 
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relativistic treatment for imperfect fluids (see also @). 
The novelty of such an approach is that particle produc- 
tion, at the expense of the gravitational held, is an irre- 
versible process constrained by the usual requirements of 
nonequilibrium thermodynamics. This irreversible pro- 
cess is described by a negative pressure term in the stress 
tensor whose form is constrained by the second law of 
thermodynamics 1 . This macroscopic approach has also 
microscopically been justified by Zimdahl and collabora- 
tors through a relativistic kinetic theoretical formulation 
[24) . In comparison to the standard equilibrium equa- 
tions, the irreversible creation process is described by 
two new ingredients: a balance equation for the parti- 
cle number density and a negative pressure term in the 
stress tensor. Such quantities are related to each other in 
a very definite way by the second law of thermodynam- 
ics. Since the middle of the nineties, several interesting 
features of cosmologies with creation of cold dark matter 
and radiation have been investigated by many authors 

[HI SI [M III, 11 . 

In what follows, as theoretically predicted by inflation 
and observationally indicated the angular power spec- 
trum of the temperature fluctuations, we consider the 
EFE for a homogeneous, isotropic, spatially flat universe 
with gravitationally induced particle production: 



8ttG P = 3- 



8ttG( P + Pc ) = -2- - 

a a z 



(4) 



(5) 



where an overdot means time derivative, p and p are the 
dominant energy density and pressure of the cosmic fluid, 
respectively, and p c is a dynamic pressure which depends 
on the particle production rate. Special attention has 
been paid to the simpler process termed "adiabatic" par- 
ticle production. It means that particles and entropy are 
produced in the space-time, but the specific entropy (per 
particle), a — S/N, remains constant [8(. In this case, 
the creation pressure reads 04HI 



Pc 



(p + p)r 

3H ' 



(0) 



where T with dimensions of (rime) -1 is the particle pro- 
duction rate and H = a/a is the Hubble parameter [13J. 
In principle, the quantity T should be determined from 
quantum field theory in curved spacetimes by taking into 
account that particle production is an irreversible pro- 
cess. 

But how the evolution of a(t) affected by T? By as- 
suming a dominant cosmic fluid satisfying the equation 



of state (EoS), p — top, where lu is a constant, the EFE 
imply that 



H+-{l + u)H 2 (l-—\=Q. 



(7) 



The de Sitter solution (H = 0, F = 3H = constant) is 
now possible regardless of the EoS defining the cosmic 
fluid. Since the Universe is evolving, such a solution is 
unstable, and, as long as V <C 3H, conventional solu- 
tions without particle production are recovered. From 
the above equation, one may conclude that the main ef- 
fect of r is to provoke a dynamic instability in the space- 
time thereby allowing a transition from a de Sitter regime 
(r ~ 3H) to a conventional solution, and vice versa (see 
subsections A and B below). 



A. From an early de Sitter stage to the standard 
radiation phase 

Let us first discuss the transition from an initial de 
Sitter stage to the standard radiation phase. The main 
theoretical constraints are: 

• The model must not only solve the horizon problem 
but also provide a quasiclassical boundary condition 
to quantum cosmology (a hint on how to solve the 
initial singularity problem). 

• Massless particles cannot be quantum-mechanically 
produced in the conventional radiation phase 
(Parker's Theorem). 

To begin with, let us assume a radiation dominated Uni- 
verse (u> = 1/3, r = r r ). The dynamics is determined by 
the ratio T r /3H [see Eq. The most natural choice 

would be a ratio which favors no epoch in the evolution of 
the Universe (T r /3H — constant). However, the parti- 
cle production must be strongly suppressed, T r /3H <C 1, 
when the Universe enters the radiation phase. The sim- 
plest formula satisfying such a criterion is linear, namely: 
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3H 



H 



(8) 



where Hj is the inflationary expansion rate associated 
to the initial de Sitter (H < Hj). It is worth to notice 
that such particle creation rates have been previously dis- 
cussed by several authors (see @, H, H3, HI] and references 
therein). It is also worth noticing that for "adiabatic" 
photon creation the form of the blackbody spectrum is 
preserved in the course of the expansion [291 ] . 
Now, inserting Eq.® into Eq. ([7]) it becomes: 



1 The quantum approach is unable to provide the entropy burst 
accompanying the particle production since it is adiabatic and 
reversible. 
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The solution of the above equation can be written as 
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H(a) 



1 + Da 2 ' 



(10) 



where D > is an integration constant. Note that 
H = Hi is a special solution of Eq. ((9]) describing the 
exponentially expanding dc Sitter space-time. This solu- 
tion is unstable with respect to the critical value D = 0. 
For D > 0, the universe starts without a singularity and 
evolves continuously towards a radiation stage, a ~ t 1 / 2 , 
when Z?a 2 3> 1. By integrating Eq. (fit))) , we obtain the 
scale factor: 



a A / a 
fljt = In — + — — 

a* 2 V a* 



(11) 



where A 2 = Da 2 is an integration constant and a* defines 
the transition from the de Sitter stage to the beginning 
of the standard radiation epoch 2 . At early times (a <C 
a*), when the logarithmic term dominates, one finds a ~ 
a r e H ' f , while at late times, a ^> a^,H <C -ff/, Eq. (fTTj) 

1 /2 

reduces to a ~ a* (^rs^t) > and the standard radiation 
phase is reached. 

It should be noticed that the time scale Hi^ 1 provides 
the greatest value of the energy density, pi = 3H] / 8nG, 
characterizing the initial de Sitter stage which is sup- 
ported by the maximal radiation production rate, r r = 
3Hi. From Eqs. (j3| and (fTU|) we obtain the radiation 
energy density: 



Pr = Pi 



1 + A 2 



(12) 



As expected, we see again that the conventional radiation 
phase, p r ~ a~ , is attained when a ^> a*. 

Now we pose the following question: how does the cos- 
mic temperature evolve at the very early stages? For 
"adiabatic" production of relativistic particles the energy 
density scales as p r ~ T 4 [1, [t| , and the above equation 
implies that 



see that the expansion proceeds isothermally during the 
de Sitter phase (a <C a*) which means that the super- 
cooling and subsequent reheating that is taking place in 
several inflationary variants are avoided [34|,[35|. In other 
words, there is no 'graceful exit' problem. 

After the de Sitter stage, the temperature decreases 
continuously in the course of the expansion. For a 3> 
a* (H <C Hi), we obtain T ~ a -1 . Accordingly, the 
comoving number of photons becomes constant since n oc 
a -3 , as expected for the standard radiation stage. 3 

In this context, we also need to answer to the following 
question: what about the initial temperature Tjt Since 
the model starts as a de Sitter space-time, the most nat- 
ural choice is to define Tj as the Gibbons-Hawking tem- 
perature [36| of its event horizon, Tj = Hi /2ir. Naively, 
one may expect Tj of the same order or smaller than 
the Planck temperature because of the classical descrip- 
tion. From EFE we have pi — 3mp 2 Hi 2 /8tt (where 
mpi ~ 1.22 x 10 19 GeU), and since the energy density is 
pi = N*(T)Tf, one finds I> ~ Hi ~ 10 19 GeV (where 
N*(T) = ir 2 g*(T)/30 depends on the number of effec- 
tively massless particles). 

Naturally, due to the initial de Sitter phase, the model 
is free of particle horizons. A light pulse beginning at 
t = — oo will have traveled by the cosmic time t a physical 
distance, <fe(i) = a(t) f ^ -~, which diverges thereby 
implying the absence of particle horizons. The latter fea- 
ture means that the local interactions may homogenize 
the whole Universe. 

Since photons are not produced in the radiation phase, 
the Big-Bang Nucleosynthesis (BBN) may work in the 
conventional way [37j . Subsequently, the Universe en- 
ters the cold dark matter [Einstein-de Sitter, a(t) oc i 2 / 3 ] 
dominated phase. Finally, we have also verified that a 
large class of r r is capable to overcome the "graceful 
exit" and fine tuning problem. Indeed we have found 
that the correct transition from an early de Sitter to the 
radiation phase is valid even for T r oc H n with n > 2 
(work in progress). The details are of course different 
depending on the power n, but the qualitative fact of the 
transition is universal, and to our opinion this is very 
good news because it shows that it can be a clue for a 
general graceful exit mechanism. 



T = Ti 



1 + A 2 ( - 

a* 



-1/2 



(13) 



where Tj is the temperature of the initial de Sitter phase 
which must be uniquely determined by the scale Hi . We 



2 This kind of evolution was first discussed by G. L. Murphy I3l|l by 
studying possible effects of the second viscosity in the very early 
Universe. Later on, it was also investigated in a more general 
framework involving cosmic strings by J. D. Barrow |3^ | who 
coined the expression "Deflationary Universes" . It has also been 
discussed in connection with decaying A(t)-models [33| . 



B. From Einstein-de Sitter to a late time de Sitter 

stage 

Due to the conservation of baryon number the remain- 
ing question is the production rate of cold dark matter 
particles and the overall late time evolution. In other 
words, what is the form of r<2 m ? For simplicity, we con- 
sider here only the dominant CDM component. 



3 Since n r oc T 3 , the average photon concentration reads n r 
m[l + X 2 (a/a t ) 2 ]- 3/2 . 
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In principle, Tdm should be determined from quantum 
field theory in curved spacetimes. In the absence of a rig- 
orous treatment, we consider (phenomenologically) the 
following fact HiS]: 
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FIG. 1: Overall evolution of the scale factor predicted by the 
matter creation model (solid line) and the traditional ACDM 
cosmology (open points). In this plot we have adopted the 
best fit, fl m = 0.273, from WMAP7 [SJ. 



• All available observations are in accordance with 
the ACDM evolution both at the background and 
perturbative levels. 

Now, we recall that a flat ACDM model evolves like: 



H 



-H l 



1 - 



H 



= 0, 



(14) 



where Hf 2 = A/3 sets the Hubble scale of the final de Sit- 
ter stage (H > Hf). Such behavior should be compared 
to that predicted for a dust filled model (uj = 0, T = Tdm) 
with particle production [see Eq. 



H 



-H 2 1 



3H 



(15) 



we see that the same 
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By comparing Eqs. (|14|) and (|15 

background evolution requires that Td m /3H = (Hf / H) 
Thus, in the background solution, the particle creation 
rate does not depend on a given scale but rather it is con- 
sidered homogeneous. The limiting value of the creation 
rate, Tdm — 3i?/, leads to a late time de Sitter phase 
(H = 0, H = Hf) thereby showing that the de Sitter 



solution now becomes an attractor at late times. With 
this proviso, the solution of Eq. (|15|) reads: 

H 2 = H Q 2 E(z) 2 = H Q 2 [f2 m (l + zf + n A ] , (16) 

where il\ = (Hf/Ho) 2 = 1 — Cl m is smaller than unity 
and 1 + z = a . Such solution mimics the Hubble func- 
tion H(z) of the traditional flat A-cosmology, with Cl\ 
playing the dynamical role of f2 A (dark energy appearing 
in the concordance model) 4 . The dark matter parame- 
ter (Qdm — 1) is also replaced by an effective parameter, 
(£ld m )eff = 1 — ^a, which quantifies the amount of mat- 
ter that is clustering. This explains why this model is in 
agreement with the dynamical determinations related to 
the amount of the cold dark matter at the cluster scale, 
and, simultaneously, may also be compatible with the 
position of the first acoustic peak in the pattern of CMB 
anisotropics which requires fltotai = 1- 
By integrating Eq. (fl"6f we obtain 




sinh 3 




(17) 



Note that the late time dynamics is determined by a 
single parameter, namely: f2 A = 1 — Omi an d is identical 
to that predicted by the flat ACDM model. Using the 
current, a joint statistical analysis, involving the latest 
observational data (SNIa [4l[, BAO [H and CMB shift 
parameter |40|) is implemented. We find that the overall 
likelihood function peaks at Cl m — 0.274 ± 0.011 with 
Xt t(^m) — 543.18 for 557 degrees of freedom. Since 
the current statistical results are in excellent agreement 
with those provided by WMAP7 [lij, for the rest of the 
paper we will restrict our present analysis to the choice 
(r2 m ,CT8,o) = (0.273,0.811), where cr 8 ,o is the rms mass 
fluctuations on scales of 8/i _1 Mpc at redshift z = 0. 

In Figure 1, we show the overall evolution (radiation, 
matter and dark energy dominated eras) of our complete 
cosmological scenario which coincides exactly with the 
one recently discussed in Refs. [12|,ll5| following a slightly 
different approach. Note also that by replacing the value 
of Td m into the definition of the creation pressure [see 
Eq. ©] one obtains that it is negative and constant 
( Pc = -3H 2 f /8nG = -3n A Hg/87rG). 



III. THE EVOLUTION OF THE LINEAR 
GROWTH FACTOR 

In this section, we briefly discuss the basic equation 
which governs the behavior of the linear matter pertur- 



4 Observationally there is no reason to distinguish between Ha 
which appears in Eq. (1 1 6 I t and that of the cosmological constant 
term parametrized by Q^. 
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bations 5 m = Sp m /p m on subhorizon scales in the CCDM 
model, assuming that the particle creation rate remains 
homogeneous and only the corresponding effective dark 
matter forms structures. The reason for introducing the 
growth analysis here is to give the reader the opportunity 
to appreciate also at the perturbative level, the relative 
strength and similarities of the CCDM and ACDM mod- 
els used to constrain the growth index. As discussed by 
Jesus et al. [43| based on the Neo-Newtonian approach 
fiif , the evolution equation of the matter fluctuations of 
a CCDM cosmology reads: 



d 2 S r , 
drf 



+ F(r,)^ + g(r ] )S m = 0, (18) 



whose solution is 5 m (r]) cx T>(r]), with V{rf) denoting 
the linear growing mode (usually scaled to unity at the 
present time). The functions appearing in (|18[) are de- 
fined by: 

_, O m (l+6c p ) + 2fi A e 3 "(8 + 3c p ) 
F W = T^K , a ^ > ( 19 ) 



2(n m + tt A e 3 *>) 



9^ 



15n A e 371 {l + c p ) ~ 3» m (2 + eg) 



(20) 



where r] = lna(i) and Cl m = 1 — Oa- The quantity c p can 
be viewed as the "effective adiabatic" sound speed 



2 _ &Pc 

Ceff - ad 8 Pn 



(21) 



where 6p c is the perturbation of creation pressure. Obvi- 
ously using the above equation and Eq.([S]) we find that 
the corresponding "effective adiabatic" sound speed must 
be negative c p = <? e ** ad < 0. Note that the latter restric- 
tion is valid also for the interacting dark energy models 
(see [IE]) 5 . We would like to stress that for simplicity we 
are using a constant c p (for more discussions see (43J). In 
this context, the functions F(rj) and Q{if) defined above 
can be written in terms of Cl m (r]) as 



i + c p + 15(1 - a m ) 



(22) 



5 Based on a scalar field description (in our case see section V) 
one can show that the sound speed c 2 . = Sp^/Sp^, is exactly 
unity. Indeed using similar arguments to those of [45H we also 



prove that 
implies c 2 , 



= <5(4- + V) = . 
Valiviita et al. 



= 5(4- - V) = 8 Pcj> which 
proposed that in order to 
1 



avoid instabilities in the dark energy we have to have c 2 
and ad < 0. Obviously both conditions are full-filled here. 



and 



) = -0^ + 15(l + c p )(l-O m )-3(2 + c p )n m , (23) 



where the definition 



n,. 



n A e 3r > ' 



(24) 



has been adopted. At this point, we remind the reader 
that solving Eq. (fT8|) for the ACDM cosmology 6 , we derive 
the well-known perturbation growth factor (see (46l|): 



V A (z) 



5h m E(z) 



(1 + u)du 
E 3 {u) 



(25) 



Obviously, for E{z) ~ Vl\1 2 (1+z) 3 / 2 it gives the standard 
result T>(z) ~ a = e n = (l + z) _1 , which corresponds to 
the matter dominated epoch, as expected. 

Now, for any type of DE, an efficient parametrization 
of the matter perturbations is based on the growth rate 
of clustering originally introduced by Peebles [46( . This 



m = 



d In 8 n 
dr] 



which implies 



T> (a) = exp 



a OL{x) 



dx 



(26) 



(27) 



where 7 is the so called growth index (see Refs. (47r - 
[52[ ) which plays a key role in cosmological studies, espe- 
cially in the light of recent large redshift surveys (like 
the WiggleZ and SDSS (DR9); see and refer- 

ences therein). As an example, it was theoretically shown 
that for DE models which adhere to general relativ- 
ity the growth index 7 is well approximated in terms 
of the equation of state parameter ^de — 3 qZ-5 > ( see 
HHH El HI), which boils down to 7A « 6/11 for the 
ACDM cosmology w(z) = — 1. Notice, that in the case 
of the braneworld model of Dvali, Gabadadze <fc Porrati 
[56} we have jdgp ~ 11/16 (see also [Umi), while for 
the f(R) gravity models we have j R ~ 0.41 - 0.43 HU 
at the present time. 

Differentiating Eq. (1261) with respect to 77 we have 



(28) 



df_ + 2 _ 1 d 2 S m 
dr] ' S m dr] 2 



Using Eqs. (fTSI) . ([2^)1 and ([2^1) . we find after some algebra 
df dfl r 



dfl m dr] 



+ f + F( v )f + g{n) = 



(29) 



For the usual ACDM cosmological model we have ^d^" ~'~ 

/ 2 + X(n)f + G(ri) = where X{r,) = \ + f [1 - flZlv)} and 
g{rj) = -ffi m (j7). 
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TABLE I: Summary of the observed growth rate and references. 



Ind6x 


Z 


growth. rs,t6 (-Ao&s) 




^iittiIiaIc fin Pin- ^ 

oymuois i m r ig. o I 


i 


0.17 


0.510 ± 0.060 


[59, 60J 


open circles 


2 


0.35 


0.440 ± 0.050 


[59, 61] 


open circles 


3 


0.77 


0.490 ± 0.180 


[59, 62] 


open circles 


4 


0.25 


0.351 ± 0.058 


[53] 


open triangles 


5 


0.37 


0.460 ± 0.038 


[53] 


open triangles 


6 


0.22 


0.420 ± 0.070 


m 


solid circles 


7 


0.41 


0.450 ± 0.040 


M 


solid circles 


8 


0.60 


0.430 ± 0.040 


M 


solid circles 


9 


0.78 


0.380 ± 0.040 


M 


solid circles 


10 


0.57 


0.427 ± 0.066 


[55] 


solid square 



where 



dr) 



= -3fi m (77) [l-flm(»?)] 



Inserting the ansatz / ~ ftm^ 1 into Eq. (HHJ), using 
simultaneously Eqs. (|22|) . (|23)l . (l24l) and performing a 
first order Taylor expansion around Q m = 1 (for a similar 
analysis see [5l|, HH) we find that the asymptotic value 
of the growth index to the lowest order is 



(30) 



7 



3(13 + 12c p ) _ 7A (13 + 12cp) 
ll + 6c p 2(1 + 7ac p ) 

Inverting the above equation we have 

13 7 A - 2 7 
Cp ^2 7A ( 7 -6) ■ 



(31) 



(32) 



4 - 





FIG. 2: Le/t Pand: The variance A X 2 = X 2 



X.Tl 



around 



the best fit 7 value for the CCDM cosmology. Note that 
the cross corresponds to ( 7 a, Axi CT ) = {tt>1)- Right Panel: 
The Ax 2 versus c p . The corresponding cross is (c p , Axio-) = 
(-1.008,1). 



We have checked for various values of c p and f2 m , that 
using Eq. ([3T|) in Eq. (|2"T1) the latter provides an ex- 
cellent approximation to the numerically obtained form 
of D(t]) appears in Eq. (fT8|) . Indeed the difference be- 
tween the two approaches is less than 0.1 — 0.2%. Finally, 
from the above analysis it becomes clear that the possible 
difference between the CCDM and ACDM predictions is 
quantified only at the perturbative level, via the value of 
7, because the two cosmological models share the same 
Hubble parameter as well as the same number of free 
parameters, namely the dimensionless matter density at 
the present epoch f2 m and the growth index. In the case 
of 7 = 7a ~ 6/11 we find c* ~ —1.008. Hence expanding 
Eq. (l3Tj) around c* we can write 

6 

7 ~ IT + 7c ^ Cp ~ c ^ ' ^ 

where 

dl , ^ x 162 

7 « = ^ = WTk? ■ (34) 



IV. FITTING THE CCDM GROWTH INDEX 
TO THE DATA 

In the following we briefly present some details of the 
statistical method and on the observational sample that 
we adopt in order to constrain either the growth index 
or the "effective adiabatic" sound speed, presented in the 
previous section. 

A. The growth data 

The growth data that we will use in this work based 
on 2dF, VVDS, SDSS and WiggleZ galaxy surveys, for 
which their combination parameter of the growth rate of 
structure, f(z), and the redshift-dependent rms fluctua- 
tions of the linear density field, (J%{z) 1 is available as a 
function of redshift, f(z)a%(z). The fas = A estimator 
is almost a model-independent way of expressing the ob- 
served growth history of the universe [59| . In particular 
the data used are based on 

• The 2dF (Percival et al. [13), SDSS-LRG 
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(Tegmark et al. [611) and VVDS (Guzzo ct al. j62|) 
bas ed g rowth results as collected by Song & Perci- 
val |59( | . This sample contains 3 entries. 

• The SDSS (DR7) results (2 entries) of Samushia 
et al. [53| based on spectroscopic data of ^106000 
Luminous Red Galaxies (LRGs) in the redshift bin 
0.16 <z< 0.44. 



• The WiggleZ results of Blake et al. [54( based on 
spectroscopic data of ^152000 galaxies in the red- 
shift bin 0.1 < z < 0.9. This dataset contains 4 
entries. 

• The SDSS (DR9) results of Reid et al. [U based 
on spectroscopic data of ^264000 galaxies in the 
redshift bin 0.43 < z < 0.70. This dataset includes 
1 entry. 

In Table I we list the precise numerical values of the data 
points with the corresponding errors bars. 



B. Observational constraints 

In order to constrain the CCDM growth index (or 
Cp) we perform a standard \ 2 minimization procedure 
between the growth data measurements (see previous 
section), A obs = f bs(z)a s ^ obs (z), with the growth val- 
ues predicted by the CCDM model at the correspond- 
ing redshifts, A(p,z) = f(p,z)ag(p,z) with er 8 (p, z) = 
as,oD(p, z). The vector p contains the free parameters 
of the cosmological model. In particular, the essential 
free parameters entering in the theoretical 



where <Ji is the observed growth rate uncertainty. Note 
that we sample 7 <E [0.1, 1.3] in steps of 0.001. 

In the left panel of Fig. 2 we show the variation of 



Ax 2 



X 



x 2 ■ 



around the best 7 fit value for the 



v 2 

A mm 



CCDM cosmology. We find that the likelihood function 
peaks at 7 = 0.60 ± 0.072 (or c p = -1 ± 0.011; see the 
right panel of Fig. 2) with Xmin — 7 -^^ f° r ^ degrees of 
freedom. Notice, that for the physically acceptance range 
0.22 < f2 m < 0.32, we obtain either 0.50 < 7 < 0.70 or 

-1.016 < c p < -0.983 with XmiJ 8 G I - 90 ' L03 ]- Hence 
the effective sound speed varies very little in function 
of Cl m . Alternatively, considering the ACDM theoreti- 
cal value of 7 (= 6/11) and minimizing with respect to 
n m0 we find f2 TO = 0.243 ± 0.034 (see also [52]) with 
/dof ~ 7.37/8. 
Our best-fit 7 value is in agreement within la (Axf^ — 
1) uncertainty with the theoretically predicted value of 
7A ~ 6/11 (see the cross in Fig. 2). Finally, in Fig. 3, 
we plot the measured A b s (z) with the estimated growth 
rate function, A(z) = f(z)a$(z) (see solid line). 

We would like to finish this section with the follow- 
ing observation: if the particles are created proportional 
to the DM density (see, for instance @,H3,|2^|) then at 
each point the growth of matter density perturbations 
will speed up with respect to an homogeneous parti- 
cle creation rate, changing the growth factor and thus 
it could potentially affect the observational constraints. 
However, in our case we do not face such a problem 
because the corresponding particle creation rate term 
(^dmSH = (Hf/H) 2 ) in the matter dominated era be- 
comes Ydm/SH oc 1/p. Notice that in order to obtain 
the latter relation we use Eq.([3]). 



0.6 



0.2 



0.2 



0.4 



0.6 



V. SCALAR FIELD DESCRIPTION 

Matter creation models constitute a possible way to 
explain the cosmic acceleration without the introduction 
of a dark energy component. However, it is sometimes 
desirable to represent the cosmic evolution in a field the- 
oretical language, i.e., in terms of the dynamics of an 
ordinary scalar field (<j>). In a point of fact, all the dy- 
namical stages discussed here can be described through a 
simple scalar field model (for a similar analysis see [28] ) . 

To begin with, let us replace p and p to t = P+Pc in Eqs. 
(HI) and §5§ by corresponding scalar field expressions 



FIG. 3: Comparison of the observed and theoretical evolu- 
tion (see solid line) of the growth rate of clustering A(z) = 
f(z)ag(z). The thin-line error bars correspond to la 7 un- 
certainties. The different growth data sets are represented by 
different symbols (see Table I for definitions). 

expectation are: p = (7, fl m ). The \ 2 function is defined 
as: 



v{<t>), 



Ptot -> P</> 



V{<t>) ■ (36) 



Inserting the latter into the Friedmann's equations we 
can separate the scalar field contributions and express 
them in terms of H and H, i.e., 



1 



4ttG 



(37) 



N r 



! Np) = E 



A bs{zi) - A(p,Zi) 



(35) 



8vrG I 



H 
3IP 



3ff 2 
8^G 



aH 
~3H 



(38) 
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where H = aHH and prime here denotes derivative with 
respect to the scale factor. Now, considering that dt = 
da/aH, Eq. (|37l) can be integrated to give 



H 

4^G 



1/2 



dt 




A. Early de Sitter - radiation: Deflationary stage 



In this case the Hubble function is given by Eq. (|T0|) . 
Obviously we can integrate Eq. (|39l) . in the interval [0, a] 
to obtain 

^( a ) = sinh" 1 (y/Da 

V 2ttG v 

= = In ( y/Da + VDa 2 + . (40) 

Note that at the time of inflation (a = a*) the corre- 
sponding scalar field is 



*.= * i„(a + v^TT) 



(41) 



Also after some simple algrebra, the potential energy be- 
comes 



V{c 



Hj 3 + Da 2 



&ttG (l + Da 2 ) 2 ' 



V{4>) 



Hj 3 + sinh 2 (V2^G 



8ttG [l + sinh 2 (V2^G0)] 2 



(42) 



(43) 



In the context of slow-roll approximation, one can prove 
that the density fluctuations are of the form Sh ~ 
H 2 /cb 2 ~ (see Peacock, Cosmological Physics 1999). 
Using Eqs. (fT0|) and (|37|) the function H 2 /<fi 2 becomes 

£ = -4^£ = -4.G* = 27rG ii±^ . (44) 
02 # ai? Da 2 

At the epoch of inflation (a = a*) we get 



H\ , _ ^ (1 + A 2 ) 



(45) 



Inserting the latter into the density fluctuation constrain 
we obtain 



1 S H 1 #h 



A 2 2ttG 2tt Z? 



(46) 



(47) 



where Z p ; ~ yG is the Planck length in units h = c — 1 . 



B. Accelerating dark matter stage 

In this case the late evolution of the Hubble function is 
given by Eq. (jTB"]) . Now, inserting Eq. (|TB)) as well as its 
derivative (H ) into Eq. (1591 and integrating we obtain: 



A In 



1 - fi A + A a 3 + Vi - 0/ 



(48) 



where A = 1/y 24ttG(1 - fi A ). From Eqs. ([T5J and 
it is possible to show that 



V(a) 



Hi 



~ 1 + 2» A _ 3 

S2a a 



(49) 



BttG 

Finally, combining Eqs. (|48p and (|4"9")l we find 

7(0) = B + Gcosh(w0) , 
where B = Hf(5fi A + ^a)/ 3271 "^ 1 ~ ^a), C = "^(i 



(50) 



20\)r2 A /327rG(l - Oa) and w = a/247tG(1 - fi A ) 



VI. FINAL REMARKS 

A new cosmology based on the production of massless 
particles (in the early de Sitter phase) and CDM parti- 
cles (in the transition to a late time de Sitter stage) has 
been discussed. The same mechanism avoids the initial 
singularity, particle horizon and the late time coincidence 
problem of the ACDM model has been phenomenologi- 
cally eliminated (A = 0) because there is no dark energy 
in our accelerating scenario. In particular, this means 
that the dark energy component may be only a "mirage" 
(an effective description) , since it can be mimicked (glob- 
ally and locally) by the gravitationally induced particle 
production mechanism acting in the evolving Universe. 

In this scenario, the standard cosmic phases - a radia- 
tion era followed by an Einstcin-de Sitter evolution driven 
by nonrelativistic matter until redshifts of the order of a 
few - are not modified. However, the model has two ex- 
treme accelerating phases (very early and late time de 
Sitter phases) powered by the same mechanism (particle 
creation). Therefore, it sheds some light on a possible 
connection among the different accelerating stages of the 
universe. In particular our model, since H 2 = (1\Hq, 

where fi A ~ 0.7 and H ~ 1.5 x 10~ 42 GeU, it sets the 
ratio of the primeval and late time de Sitter scales to 
be pi/pf = {HijHj) 2 ~ 10 122 . This large number is 
ultimately obtained through a combination of thermo- 
dynamics and quantum field theory in curved spacetimes 
in virtue to the association of the Hawking-Gibbons tem- 
perature (T = H/2-k) to the early de Sitter phase. Due 
to the maximal radiation production rate, T r ^ max = SHi, 
the model deflates from an unstable de Sitter with initial 
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expansion rate (H = Hi). It is exactly this Gibbons- 
Hawking connection (and the present value of Hq) that 
lead for the given number under the proviso that the Uni- 
verse starts with the Planck temperature. In principle, 
such a result in the present context has no correlation 
with the so-called cosmological constant problem (in this 
connection see 33]). In a forthcoming communication, 
it will be shown that the above ratio is preserved even 
if the initial de Sitter phase is powered by a production 
rate term proportional to H n , n > 2. 

As it appears, the cosmic history discussed here is 
semiclassically complete. However, there is no guaran- 
tee that the initial de Sitter configuration is not only the 
boundary condition of a true quantum gravitational ef- 
fect. In other words, the very early de Sitter phase may 
be the result of a quantum fluctuation which is further 
semi-classically supported by the creation of massless 
particles (in this connection see [g3[ and Refs. therein). 

Naturally, the existence of an early isothermal de Sit- 
ter phase suggests that thermal fluctuations (within the 
de Sitter event horizon) may be the causal origin of the 
primeval seeds that will form the galaxies. Such a pos- 
sibility and its consequences for the structure formation 
problem deserves a closer investigation and is clearly out 
of the scope of the present paper. 

We stress that our model provides a natural solution to 
the horizon problem and finally it connects smoothly the 
radiation and the matter dominated eras, respectively. 
At late times it also mimics perfectly the cosmic expan- 
sion history of the concordance ACDM model. In this 



context, we discuss the behavior of the linear matter per- 
turbations 6 m on subhorizon scales for the CCDM model 
and the main results are shown in Figs. 2 and 3. 

For completeness, we have also represented the evolu- 
tion of our model in terms of the dynamics of an ordinary 
scalar field (4>) and derived analytically the scalar field 
potential for two regimes: (i) when the Universe evolves 
from an early de Sitter to a Radiation phase and (ii) when 
it goes to the CDM phase to a late time de Sitter stage. 

At present, we also know that a more complete version 
of the late time evolution must be filled with CDM (~ 
96%) and baryons (~ 4%), and, unlike A cosmology, the 
baryon to dark matter ratio is a redshift function (T3 , [l5j . 
In particular, this means that studies involving the gas 
mass fraction may provide a crucial test of our scenario, 
potentially, modifying our present view of the dark sec- 
tor. Some investigations along the above discussed lines 
are in progress and will be published elsewhere. 
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